Three-photon bosonic coalescence in an integrated tritter 
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The main features of quantum mechanics reside into 
interference deriving from the superposition of differ- 
ent quantum objects. While current quantum optical 
technology enables two-photon interference both in bulk 
and integrated systems, simultaneous interference of 
more than two particles, leading to richer quantum 
phenomena, is still a challenging task. Here we report the 
experimental observation of three-photon interference 
in an integrated three-port directional coupler realized 
by ultrafast laser-writing. By exploiting the capability of 
this technique to produce three-dimensional structures, 
we realized and tested in the quantum regime a three- 
photon beam splitter, namely a tritter, which allowed us 
to observe bosonic coalescence of three photons. These 
results open new important perspectives in many areas 
of quantum information, such as fundamental tests of 
quantum mechanics with increasing number of photons, 
quantum state engineering, quantum sensing and quan- 
tum simulation. 

The Hong-Ou-Mandel effect is a purely quantum phe- 
nomenon deriving from the interference of two photons JJ]. 
It is a consequence of the bosonic nature of light and can be 
observed when two identical photons impinge simultaneously 
on a two-port balanced beam splitter and, as a result, both 
are forced to exit from the same port. This effect can be ex- 
ploited to quantify the amount of indistinguishability of cor- 
related photons yl IH or to project a quantum state onto the 
maximally-entangled singlet Bell state, representing an im- 
portant resource in many research areas of quantum informa- 
tion, from quantum computation to quantum metrology J3-01 ■ 
The interaction of single-photon states on a beam splitter and 
successive post-selection may indeed lead to the generation of 
entangled states [0] such as NOON and GHZ states Hfl. As a 
matter of fact, achieving interference of more than two par- 
ticles on a beam splitter is a challenging task JjJ, however it 
could open the door to the observation of quantum phenomena 
in increasing size systems and to the engineering of quantum 
states for quantum information processing fioifllll . 

The adoption of a multiport beam splitter has been proposed 
to perform additional classes of possible EPR-Bell experi- 
ments in a Hilbert space with dimension higher than two 11211 . 
Furthermore the capability of implementing multipath inter- 
ferometers by multiport beam splitters could be adopted to 
perform quantum computation tasks f 14]. Recently a general- 
ized Alport beam splitter has been theoretically addressed in 




FIG. 1 : (a) Optical scheme for a three-modes splitter obtained by a 
set of conventional beam-splitters and phase shifters fl2l[l3ll . (b) In- 
tegrated 3-dimensional scheme for a tritter. (c) Three-photon bosonic 
coalescence effect in a tritter: all the terms with two photons exiting 
on the same output port vanish due to bosonic coalescence. 
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and the behaviour of N photons impinging onto 
the Alport beam splitter has been studied by using a statistical 
approach, since the complexity of this problem grows faster 
than exponentially with growing N. Very recently, the appli- 
cation to quantum interferometry of three and four port been 
splitters have been theoretically investigated, and the quan- 
tum enhancement in terms of sensitivity has been predicted 
1 1711 . However, three particle interferometry is mostly unex- 
plored territory, both experimentally and theoretically. The 
three beam extension of a beam splitter, the so called 'tritter', 
was first addressed in Ref. 11211 and its proposed implemen- 
tation consisted in the combination of two port beam splitters 
and phase shifters [see Fig. Q](a)] IU2LI13I1 . Three dimensional 
implementations of three- and four-port beam splitters, based 
on fused optical fibers Hill or femtosecond laser waveguide 
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writing 119142111 . have been proposed. However, only classical 
light or two-photon interference experiments were carried out 
on such devices. Therefore, an experimental investigation of 
quantum interference of more than two-photons is still lack- 
ing. 

Here we report on the first experimental investigation of 
three-photon interference in a tritter. The tritter is realized 
by femtosecond laser waveguide writing and it consists of 
three waveguides approaching three-dimensionally, thus en- 
abling simultaneous interaction of the three photons without 
having to decompose the process into cascaded two-mode in- 
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teractions and phase shifters [Fig. [U(b)]. The three photons, 
each coupled to a different waveguide, mutually interfere in 
the three-arm directional coupler by evanescent field interac- 
tion. The behaviour of a three-photon state entering the trit- 
ter is investigated by comparing the output state probabilities 
with the expected classical ones, evidencing the purely quan- 
tum effects involved in bosonic coalescence. 
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FIG. 2: Hong-Ou-Mandel effect for a |1, 1, 1) input state in an ideal 
tritter. (a-c) Output probabilities Pi,i,i (a), -p2,i,o (b), P3,o,o (c) as a 
function of the delays Xi — CTi of the impinging photons, (d) Output 
probabilities P^- k for three indistinguishable input photons, Pij : k 
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for three distinguishable input photons, and P i k for two indistin 
guishable and one distinguishable input photons. 



RESULTS 

Theoretical description of the tritter output probabili- 
ties. - Let us first address the action of a tritter on a generic 
input state |^>). It can be expressed by a unitary matrix U, 
mapping the input field operators a\ to the output field opera- 
tors b\ with b\ = J2i tj Una], EHO (see Supplementary 
Material). In the case of an ideally symmetric tritter, where a 
photon entering in one input port has the same probability of 
exiting through anyone of the output ports, the unitary matrix 
U reads: 




(1) 



Let us consider what happens when an ideal tritter is injected 
by a three-photon state |1, 1, 1). The corresponding output 
state after the evolution induced by U 1 reads: 



l,l,l)^ci, ll i|l,l > l)+c { 3,o I o}|{3,0 J 0}>. 



(2) 



Here Ci^i = — 1/\/3 and C{ 3 , ,o} = \/2/3, where k}) 
is the symmetric superposition of the three-photon states cor- 
responding to (i,j,k) photons exiting in the three output 
ports. We first observe that, as a consequence of bosonic co- 
alescence, all the terms of two photons exiting from the same 
output port disappear [Figs. [T](c)]. In order to fully describe 
three-photon interference we need to introduce two relative 
delays between the three bosons. We theoretically calculated 
the probability of the different output states [(1,1,1), (2, 1, 0), 
(3, 0, 0)] as a function of the relative delays 
(to = 1,2) between the three photons. These probabilities are 
represented in Figs. |2](a-c) as surfaces in a three-dimensional 
space and represent an extension of the well-known Hong-Ou- 
Mandel curves. This representation enables a deeper insight 
into three-photon interference, evidencing non-trivial features 
as a function of the two delays. In particular, we observe the 
presence of three regions: (i) the three input photons are indis- 
tinguishable {Pfj k ) when xi ~ X2 ~ 0, this leads to a three- 
photon coalescence effect, resulting in a dip or a peak depend- 
ing on the detected output state; (ii) only two of the three input 

(2) 

photons are indistinguishable {P^j k ) when x m ~ ^ x n or 
x m ~ x n ; (iii) the three input photons are distinguishable 
(Pi\j fc) when xi ^ x-2 ^ 0, leading to classically correlated 
outcomes. Fig. 12(d) compares the different probabilities cor- 
responding to cases (i-iii) for each output state. 

Tritter fabrication by ultrafast laser writing. - The trit- 
ter device was realized by adopting femtosecond laser writ- 
ing technology 122142411 . This technique exploits nonlinear ab- 
sorption of focused femtosecond laser pulses that determines 
permanent and localized increase of the refractive index in 
transparent materials. Waveguides are directly fabricated in 
the material bulk by translating the sample at constant speed 
along the desired path. This technique has the unique capa- 
bility of producing integrated devices with three-dimensional 
geometries. Ultrafast laser-writing technique has rapidly be- 
come a powerful tool for demonstrating new quantum inte- 
grated devices, able to perform quantum logic operations II25I1 
as well as two-photon quantum walks |2y, 27]. The fem- 
tosecond laser written tritter consists of three waveguides ap- 
proaching and interacting in the coupling region [Fig. [U(b)]. 
Thanks to this configuration equal coupling coefficients may 
be obtained, in a way that a single photon entering in any in- 
put port has the same probability to exit from any of the three 
output ports. This three-dimensional approach, enabled by 
femtosecond laser writing, avoids the use of cascaded beam 
splitters and phases which are difficult to control and provides 
a compact device that can be effectively scaled in more com- 
plex architectures. Although conceptually simple, the tritter 
needs accurate dimensioning. As a first order arrangement 
we consider the waveguides lying at the vertex of an equilat- 
eral triangle, where key parameters are waveguide spacing and 
length in the interaction region. Actually these two parame- 
ters are not independent since, in order to achieve and equal 
output distribution, the coupling coefficient k and interaction 
length L needs to satisfy the relation kL = 2tt/9 (see Sup- 
plementary Material). Therefore, sizing of the tritter requires 
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FIG. 3: Two-photon experimental reconstruction of the unitary matrix U r of the implemented tritter compared with the theoretical one U l '. 
(a)-(c) Experimental two-photon Hong-Ou-Mandel interference for input states: (e) |1, 1, 0), (f) |1, 0, 1), (g) |0, 1, 1). For each input state all 
possible output ones are acquired, (d)-(e) Real and imaginary parts of U r . (f)-(g) Real and imaginary parts of U . 



choosing the waveguide spacing d to the smallest value that 
avoids overlap (this will keep the device as short as possible). 
The chosen value of d in turn defines the coupling coefficient 
k and thus, according to the previous relation, the interaction 
length L. However, the slightly elliptical guided mode of the 
femtosecond laser written optical waveguides and other fab- 
rication imperfections require performing a second order fine 
tuning of the device that can be accomplished by slightly shift- 
ing one of the three waveguides in a scalene triangle geometry. 
This fine tuning is performed most efficiently by an iterative 
procedure where we take full advantage of the rapid prototyp- 
ing capabilities of femtosecond laser waveguide writing (see 
Supplementary Material for the final optimized geometry). 

Two-photon characterization of the tritter device. - To 
characterize the implemented device and to compare its uni- 
tary evolution with respect to the ideal case, we implemented 
the reconstruction method proposed in Ref . ll28ll . This tech- 
nique consists in measuring all the possible two-photon Hong- 
Ou-Mandel interference fringes and reconstructing the ele- 
ments of the evolution matrix U r through a minimization pro- 
cedure (see Supplementary Material for more details). We 
show in Fig. [3] the real and imaginary parts of the recon- 
structed matrix U r at 795 nm wavelength, compared to the 
corresponding theoretical matrix U f . To quantify the qual- 
ity of the fabrication process we evaluated the similarity S 
between the two corresponding Hong-Ou-Mandel visibility 
matrices, defined as S r - t = 1 - Ylijtj J2k^i l(V*kj;M ~ 
(V)ij ;fe ,z|/18. The obtained value = 0.9768 ± 0.0004 
demonstrates a high correspondence of the implemented de- 
vice with an ideal tritter. 

Experimental three photon interference. - We experi- 
mentally tested the three-photon coalescence by adopting the 
experimental setup shown in Fig. @] Four photons are pro- 
duced by parametric down conversion: one of them acts as 
the trigger for coincidence detection, while the other three are 



coupled inside the tritter, after passing through different delay 
lines. The output modes are detected by using optical fiber 
beam splitters and single-photon avalanche photodiodes. Co- 
incidences between different detectors allow us to reconstruct 
by post selection the probability of obtaining a given output 
state. By changing the relative delays between the three pho- 
tons it is possible to observe the coalescence effects in the 
different experimental conditions, as shown in Figs. |2](a-c). 
The tritter device, first characterized by single and two-photon 
states at wavelength 795 nm (see Fig. [3]) has then been tested 
by three photons at wavelength 785 nm, due to availability 
of higher transmissivity spectral filters at that wavelength. It 
has to be noted that the unitary matrix of the tritter W at 785 
nm is still in very good agreement with that of an ideal trit- 
ter (see Eq. [B, yielding a similarity S = 0.9595 ± 0.0003 
(see Supplementary Material). We have injected a three pho- 
ton state |1, 1, 1) and observed the different output contribu- 
tions |1, 1, 1), |2, 1,0), |3, 0,0). Due to the distinguishability 
between photons belonging to different pairs the expected vis- 
ibilities of peaks (or dips depending on the observed contri- 
butions) change depending on which photon is delayed. We 
compare in Fig. [5] the experimental and theoretical output 
states probabilities, obtained for different delayed photons at 
fixed output state. The expected probability values have been 
calculated by introducing in the theoretical model the distin- 
guishability parameter p, which quantifies the distinguisha- 
bility between the two different pairs, and by including six- 
photon events generated by the source. We investigated three- 
photon interference when delaying one of the two indistin- 
guishable photons [case (A) in the figure] or when delaying 
the partially distinguishable one [case (B) in the figure]. In 
addition, we also considered the case when the partially dis- 
tinguishable photon is kept out of the interference region and 
the two indistinguishable ones have a relative delay [case (C) 
in the figure]. The latter case corresponds to two-photon inter- 
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FIG. 4: Experimental setup for the characterization of the tritter and the three-photon coalescence experiment. Two-photon and three-photon 
states, generated by parametric down-conversion, are injected in the tritter device after the propagation through spatial delay lines. Then, 
coincidence detection at the output ports are performed to reconstruct the probability of obtaining a given output state realization (i, j, k). 



ference although three photons are input and measured. It is 
observed that cases (A) and (B) are clearly different from case 
(C). This demonstrates that the output states that we observe in 
cases (A) and (B) are indeed due to three-photon interference 
and cannot be attributed to two-photon interaction. We also 
compared the obtained visibilities with the classical bounds, 
calculated by injecting the tritter with three equal-amplitude 
coherent states with randomized phases ll29ll . We observe that 
the results obtained with the |1, 1, 1) input state outperform 
classical predictions, thus demonstrating the quantum origin 
of the measured three-photon interference. Details on the data 
analysis can be found in the Supplementary Material. 
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DISCUSSION 

We reported on the experimental observation of a three- 
photon bosonic coalescence effect occurring within an inte- 
grated, 3-dimensional, 3-port beam splitter realized by the 
ultrafast laser-writing technique. To the best of our knowl- 
edge, this result represents the first experimental observation 
of the bosonic coalescence of three undistinguishable photons 
in a symmetric three ports device. The intrinsically stable and 
compact nature of this device makes it a fundamental build- 
ing block for future complex networks implemented by inte- 
grated optical circuits. These devices are expected to open 
new perspectives in many research areas of quantum informa- 
tion. Indeed, the Alport beam splitters may find a wide range 
of applications in both quantum interferometry and quantum 
metrology. A first theoretical analysis performed in Ref. Il7ll . 
shows that these devices can be exploited to build the A-mode 
interferometers where the adoption of quantum input states 
can lead to a significant enhancement in phase estimation pro- 
tocols. Furthermore, the capability of conditionally generat- 
ing path-entangled states may be increased by adopting struc- 
tures of A-port beam splitters rather than conventional two- 



FIG. 5: Visibilities V = (Too - To) /Toe for different choices of 
the delayed photon and of the measured output state contribution. 
Here, Too and To correspond respectively to the number of four-fold 
coincidences when the delayed photon is respectively out or inside 
the interference region. (A) Identical photon on input port 2 delayed 
(three-photon interference), (B) distinguishable photon on input port 
3 delayed (three-photon interference), (C) distinguishable photon on 
input port 2 out of interference region and identical photon on in- 
put port 3 delayed (two-photon interference). Blue bars: theoretical 
prediction obtained from the reconstructed tritter matrix W and par- 
tially distinguishable photon-pairs. Red bars: experimental results. 
Cyan bars: expected visibilities with classical states. 



mode couplers 13011 . This is expected to lead to a correspond- 
ing reduction in complexity of the experimental schemes in 
terms of number of optical elements. 

Several other contexts may benefit from the adoption of Al- 
port beam splitters. For instance, they can be adopted for the 
realization of "proof-of-principle" quantum simulators 13111 . 
and, combined in a modular structure, they can be used to 
realize full-scale quantum simulators for large size quantum 
systems. The adoption of multiphoton-multimode platforms 
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may indeed disclose the "hard-to-simulate" scenario by adopt- 
ing linear optics to implement a computational power beyond 
the one of a classical computer B2I1 . Furthermore, the devel- 
opment of complex integrated quantum devices may lead to 
implement fundamental tests of quantum mechanics, such as 
nonlocality tests, for quantum systems of increasing dimen- 
sionality 133 1. 
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Methods 

Waveguides fabrication.- Waveguides are directly written 
in borosilicate glass substrate (EAGLE2000, Corning Inc.) 
using a cavity-dumped Yb:KYW oscillator, which provides 
300 fs pulses at 1 MHz repetition rate. The laser beam is 
focused by a 0.6 NA 50X microscope objective 170 /im (av- 
erage) under the glass surface. Optimum irradiation parame- 
ters for single-mode waveguiding at 795 nm wavelength are 
220 nJ pulse energy and 40 mm/s translation speed. Sam- 
ple translation at uniform speed is provided by high preci- 
sion, three-axis air-bearing stages (Aerotech FiberGlide3D, 
Aerotech Inc). The produced tritter is subsequently pigtailed 
by bonding fiber arrays on both sides of the device in order to 
obtain a compact and portable device. 

Four-photon source.- Let us refer to the experimental 
setup shown in Fig. [4] in the main text. A 2 mm thick (3- 
barium borate crystal (BBO) cut for type-II phase matching, 
is pumped by the second harmonic of a Ti:Sa mode-locked 
laser beam, and generates via spontaneous parametric down 
conversion at second order two polarization photon pairs over 
two different spatial modes. The four-photon terms 1134 13511 
of the generated state are spectrally filtered by interferential 
filters (wavelength A = 785 nm, AA = 3 nm) and selected by 
four single mode fibers placed at the output of two polarizing 
beam splitters. The spatial and temporal walk-off is compen- 
sated by inserting a waveplate and a 1 mm thick BBO crystal 
on each spatial mode. One of the four photons is detected to 
work as trigger of the experiment. The other three photons 
are injected each into a single mode fiber and, after passing 
through three delay lines and after polarization compensation, 
are coupled inside the tritter device. The output modes are 
split in three equal parts by adopting fiber beam splitters and 
are detected by single photon avalanche photodiodes. Coin- 
cidences between different detectors allows us to reconstruct 
the probability of obtaining a given output state realization. 
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I. THEORETICAL DETAILS 



A. Mathematical description of tritter 



The three-photon interference phenomena in a balanced integrated tritter can be modelled by considering the propagation 
of three photons in three equally-coupled waveguides. If Aj are the annihilation operators for a photon in the j-th waveguide 
(J = 1,2, 3), their evolution along the propagation coordinate z, in the Heisenberg picture, is described by the following 
equations Qj: 



d4 

dz 
dz 



«4 



(i) 



= kA\ + kA\ + (3A\ 



where j3 is the propagation constant and k is the coupling coefficient between neighbouring waveguides (assumed constant for 
each couple of waveguides). 

Equations (fl]i admit an explicit solution as a function of z: 

k4(z), 

with: 
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One can easily calculate that a balanced tritter, which distributes each input photon with equal probability to the three outputs, 
is obtained for a propagation length z — L such that kL = |7r. For compactness of notation we can define aj = 4(0) as me 

creation operator of a photon at the j-th input port of the tritter and &t = (L) as the creation operator of a photon at the j-th 
output port of the tritter. It holds: 



U 



and 



(4) 



with: 




(5) 
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This matrix turns out to be equivalent, up to multiplications for a phase factor of rows and columns, to the real-bordered unitary 
matrix U*: 



1 



f l 1 1 

J/t = t I \ e »27r/3 g447r/3 
V3 I ^ e i47r/3 e 487r/3 



(6) 



B. Generalized bosonic coalescence 



By following the formalism introduced in Ref. [01 we now analyze the generalized coalescence effect in a configuration of m 
(to = 2,3) bosons prepared in the 3 input ports of a tritter. 

Let us first consider different three-photon input states, and compare the expected output probabilities relative to a quantum 
and a classical description: 
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We observe that for the |1, 1, 1) and |2, 1, 0) input states three-photon quantum interferences is responsible of both enhancement 
and suppression of the output probabilities. 

A similar result is obtained when two photons are impinging on the device. The comparison between the quantum and classical 
probabilities to obtain a given state at the exit of the three-port beam splitter is reported in the following table: 
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C. Multimode theory for three-photon bosonic coalescence 

We now consider the case of three independent indistinguishable photons, each of them described by the wavepacket operator: 

/>oo 

a\{n) = / dUifii^ati^e-^, (7) 



where i — 1, 2, 3 labels the optical modes k^, t; is the time delay due to spatial propagation, fi(uii) is the spectral function of 
the photon wavepacket with central frequency ljq and width Slj: 
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The output state after propagation through a tritter is obtained by replacing the time evolution equations a\(u)i) — J^i j(M 
in the input state: 



Xl«l(^)|o> = 
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FIG. 1. a Cross-section of an ideal three- waveguide coupler. The waveguides are disposed on the vertices of an equilateral triangle, with 
reciprocal distance d. b Cross-section of the three-waveguide coupler (tritter) realized in this work. To compensate for various non-idealities 
of the waveguide modes and fabrication process, the position of the third waveguide is tuned in the two directions. The waveguides constitute 
the vertices of a scalene triangle; the ideal equilateral geometry is slightly deformed. 



The probability Pi. 1,1 of obtaining one photon on each output port is calculated by selecting the terms 6{(u; m ) 62(^)63(0^) of 
the output state. For the symmetric tritter \T\ = \R\ = 1/V3, we obtain: 

Pl,l,l(Tl,T 2 ,T 3 ) = -[2 - e - Su2(Tl - T2)2 - e -^ 2 (ri-r 3 ) 2 _ (r 2 -r 3 ) 2 + 4e -fc 2 (r 1 2 +r 2 2 +r ; ?-r 1 r 2 -r 1 r3-r 2 r 3 )^ (1Q) 

By an analogous procedure, the probabilities i-2,1,0 an d p3,o,o of obtaining the outcomes (2, 1, 0) and (3, 0, 0) read: 
and: 

P 3 ,0,0(Tl, T 2 ,T 3 ) = -L[l +e -<^Vi-^) 2 +e -Su 2 (r 1 -T 3 f +fr fc 2 (r 2 -T 3 ) 2 +2e -5 W 2 (r 1 2 +r 2 2 +r|-r 1 r 2 -r 1 r3-T 2 r 3 )^ (12) 

We observe the presence of the two-photon interference terms e~ <5u2 ( Tm ~ r ™' 2 , depending only on the relative delay between 
input photons m and n, and a pure three-photon interfence term e~ 5uj ' Tl + r 2+ T 3 - t i t 2-t 1 t 3 -t 2 t 3 ) ^ d e p enc h n g on the relative 
delay between the three input photons. 

We conclude by discussing the effect of spectral mismatch between the photons in the Hong-Ou-Mandel output probabilities 
Pi.j,k(T~i, T2, T3). In this case, we can expect that the two-photon interference terms will be reduced by a factor S mt7l , describing 
the spectral overlap between the modes m and n. Analogously, the three-photon interference term will be reduced by a factor 



II. GEOMETRY OPTIMIZATION OF THE FEMTOSECOND LASER WRITTEN TRITTER 

Ideally, if waveguides with perfectly circular cross-sections and guided modes were available, a symmetric coupling between 
the three waveguides of a tritter could be achieved by means of a symmetric three-dimensional geometry. In other words, the 
three waveguides should stay at the vertices of an equilateral triangle (Fig. Q]a) m the coupling region (i.e. in the region in which 
they are close to each other and couple by evanescent field). In fact, the coupling coefficient for two neighbouring identical 
waveguides would only depend on the relative distance d. In the ideal case, as shown in Section ITAl once distance d is fixed it 
would be sufficient to adjust the interaction length L to obtain a balanced tritter. 

As a matter of fact, the femtosecond laser written waveguides that we use in this work support elliptical guided modes [3], 
which give an angle-dependent coupling [4], and other fabrication non-idealities may occur (see e.g. Ref. Hence, in order 
to realize an accurately balanced tritter, it is crucial to develop a robust experimental procedure for optimizing the fabrication 
parameters, which can compensate all the non-idealities. 

As a first step we chose the separation d between the two waveguides 1 and 2, which are at the same depth (see Fig. QJi). The 
coupling coefficient n exponentially increases[4] with decreasing d: a small value of d is preferable to obtain a compact device. 
In our case it was chosen d — 1 1 /im, which allows to fabricate three waveguides with the geometry of Fig. QJ without touching 
each other for every vertical distance h of the third waveguide. 
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The choice of d determines also the value of k. As shown in Section lTAl in a balanced tritter, kL = |7r: also L is fixed by the 
choice of d. Actually, the interaction length L can be optimized more conveniently by studying two-waveguide couplers, which 
are much simpler devices. For a two-waveguide coupler, equations similar to ([TJ can be written [6] and the coupler reflectivity^ 
results to be: 

R = sin 2 (kL) . (13) 

Hence, we fabricated several devices consisting of the waveguides 1 and 2 at distance d with different lengths L, but without 
the third waveguide i.e. simple planar directional couplers: optimum interaction length resulted L ~ 2.5 mm, which gives a 
reflectivity R = sin 2 (§77) ~ 41%. 

The next step is the optimization of the position of the third waveguide. First, it is important to note that the ellipticity of the 
guided modes, which gives an angular-dependent coupling, requires the distance between waveguides 1 and 2 to be different 
with respect to the distances between the same waveguides and waveguide 3. The former equilateral triangle on whose vertices 
the waveguides should stay in the ideal case, should now become isosceles. 

Actually, we experimentally observed slight left-right imbalance of the output splitting ratio of tritters fabricated with isosceles 
triangle geometry. We also observed that the writing order of the three waveguides influenced this imbalance. This means that 
the inscription process of a second waveguide in close proximity of another one, already fabricated, may perturb slightly the 
first one or, alternatively, be slightly disturbed by it. One can compensate this asymmetry by a small translation A (Fig. QJ>) 
of the third waveguides. Thus, a better balance between the coupling coefficients can be obtained by adjusting the position of 
the third waveguide not only in the vertical direction (distance h in Fig. [TJ)), but also in the horizontal direction (distance A 
in the same figure): the waveguides are on the vertices of a scalene triangle. We experimentally performed this optimization 
by fabricating several tritters (with d and L fixed from the previous optimization) with different values of h and A. The output 
power repartition was measured for coherent light injection in the different input ports. This allowed to interpolate the values of 
h and A which gave a balanced coupling. 

A final adjustment was performed again on the interaction length L, both because of day-to-day reproducibility issue of the 
fabrication process, and because the presence of the third waveguide also may perturb the other two. 

The tritter chosen to be employed in the quantum experiments yield the following geometric parameters: d = 11/im, L = 
2.9 mm, h = 10.3 fim, A = 0.9 ^m. 

Note that this optimization process required the fabrication of tens of tritters and couplers, with different geometric parameters. 
Given the fast processing speeds of femtosecond laser micromachining and since it doesn't require costly lithographic masks, 
the whole procedure resulted extremely cost-effective and required only a few days of work. 



III. RECONSTRUCTION OF THE UNITARY MATRIX OF THE TRITTER 

Here we report the reconstruction of the tritter's transition matrix. The present device has been optimized to work at A = 795 
nm, and hence we performed its characterization at this wavelength. The experimental reconstruction of the transition matrix 
of the implemented tritter has been obtained by adopting a method based on the measurement of all possible combinations of 
two-photon interference terms [8[]. 

Let U be the transfer matrix of the chip, whose elements Uij describe the transition amplitude of a photon entering in input 
port i and exiting from output port j. In the classical case the probability to detect a photon out of each output k and I, when two 
photons are injected in the inputs i and j where i ^ j = 1, 2, 3 and k 7^ I = 1, 2, 3, is given by: 

(a) (6) 

where Ui j G C is the matrix element (i, j) of U. The two terms (a) and (b) correspond respectively to: 

(a) the photon entering in i goes out of I while the photon entering in j goes out of k, 

(b) the photon entering in i goes out of k while the photon entering in j goes out of I, 

The above expression describes the classical case, since the overall probability is given by the sum of the probabilities for cases 
(a) and (b) and no interference occurs between the two paths. 

For two indistinguishable photons entering in the input ports i and j, the probability to detect a photon in the output ports k 
and I is given by: 



(15) 
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corresponding to quantum interference between paths (a) and (b). By varying the relative time delay between the two impinging 
photons, we observe the Hong-Ou-Mandel effect leading to peaks (or dips) in the output coincidence events. The visibility V of 
the non-classical peak (or dip) is given by: 

pC _ pQ 

V iJik ,i = ' J;fc p C W , (16) 
where positive values indicate a dip and negative values indicate a peak. 




FIG. 2. Visibilities of the dips at 795 nm, for the input modes and the outputs (k, I): (a) measured visibilities, (b) visibilities obtained 
with IA T and (c) visibilities expected with an ideal tritter U l . 



As shown in Ref. |8[], the two-photon Hong-Ou-Mandel effect can be exploited to reconstruct the transition matrix of the 
tritter. We model the transition matrix as: 

(\UxA |w 1|2 | |Wi, 3 | \ 

Ui,j = \U 2>1 \ \U 2 , 2 \e^ \U 2t3 \e*>>> . (17) 
\|W 3 ,i| |W 3 ,2|e^ 3 ' 2 \U 3>3 \e^J 

Here, only four free phases (from a global number of nine) are allowed to be different from zero. Indeed, a global phase can be 
neglected and other four phases can be set to zero by multiplying rows and columns for a phase factor, since the phase of the 
input state is undefined and photon-counting measurements are performed. The reconstruction then proceeds in two steps. 

Step (1). - The squared moduli of the matrix elements ■ | 2 are measured by feeding the tritter with single photons in only 
one input port (see table|I|. For each input i we measured the count rate riij in each output port j. The squared moduli of the 
transition matrix elements can be calculated as: 

H, 3 \ 2 = -^^. (18) 





Output 


Input 


1 2 3 


1 


0.37 0.34 0.29 


2 


0.33 0.32 0.37 


3 


0.28 0.38 0.35 



TABLE I. Measured single coupling coefficients 1 2 °f me tritter. 



The phases of the starting matrix are chosen to be equal to the theoretical values reported in Eq. (|6). 

Step (2). - In order to fully reconstruct the matrix U, we measure the two-photon Hong-Ou-Mandel visibilities V™ . fe ; 
corresponding to all possible combinations of two-fold correlations for all input states (|li, la, O3), |li, O2, 13) and |0j., I2, 13), 
where 1,2,3 are the input ports) [see Fig. 121(a)]. Then, we minimize numerically the distance between the theoretical visibility 
matrix V^. fe ; , calculated as a function of the transition matrix elements, and the measured visibility matrix V™ . fe ; : 

3 3 (yr . - YLuhi] 
RMS= £ £ 1 ^ 4 I . (19) 
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Here p is the factor which corresponds to the purity of the two-photon source, and a m is the error matrix associated to the 
measured visibilities. The factor p is evaluated by measuring the two-photon Hong-Ou-Mandel interference visibility in a 
symmetric 50-50 beam-splitter, which in our case corresponds to the value p = 0.94. The starting point Uq for the matrix U r is 
given by adopting the measured splitting ratios of Eq. d20b and by setting the phases e 1 ^ 1 - j to their theoretical values, leading to: 

/0.608 0.583 0.538 \ 

U 7 = 0.566 0.557e l27r / 3 0.608e l4jr / 3 . (20) 
\0.529 0.616e l4,r / 3 0.583e l87r/3 / 

The reconstructed matrix is shown in the main text and takes the following form: 

/ 0.593 ± 0.001 0.5928 ±0.0008 0.5444 ± 0.0006 \ 

U r = 0.5489 ± 0.0005 (0.5811 ± 0.0006)e I ( 2123±0 ° 2 ) (0.6008 ± 0.0006)e- l ( 2 030±0 002 ) . (21) 
\0.5886± 0.0009 (0.5575 ± 0.0007)e~ l ( 2167±0 002 ) (0.5853 ± 0.0006)e l ( 2110±0 002 ) / 

The errors on the parameters of W has been evaluted by performing M = 40 times a numerical simulation of the reconstruction 
process with random gaussian noise on the input visibilities VJ™. fc l , of variance equal to the experimental error er™ . fc t - As a 
figure of merit to quantify the quality of the algorithm we adopted the similarity S r<m — 1— Sfc^z \P^l j-k i~^Tj-k i\/(^ n ) 
between the measured visibility matrix V m [Fig. 12(a)] and the visibility matrix V [Fig. 12(b)] . We obtained a value S r,m = 
0.989 ± 0.002. The correspondence between the reconstructed matrix W and the ideal matrix U 1 of Eq. © is evaluated from 
the similarity 5 r , t = 1 - Y,kM \ V lj;k,l ~ V t j;k,l l/( 2n )' leading to: 5 r , 4 = 0.9768 ± 0.0004. These results show that the 

reconstructed matrix is in high agreement with the theoretical one. 



A. Reconstruction for input wavelength A = 785 nm 



The three-photon interference experiment described in the main text was performed by sending input states at 785 nm. In 
order to verify the quality of the tritter at this wavelength, we performed the characterization of the device at this wavelength by 
adopting the same method described in the previous section. 




FIG. 3. (a-c) Two-photon Hong-Ou-Mandel dips for input states: (a) |1, 1, 0), (b) |1, 0, 1) and (c) |0, 1, 1). (d-f) Corresponding visibilities 
of the dips at 785 nm, for the input modes (i, j) and the outputs (k, I): (d) measured visibilities, (e) visibilities obtained with U lr and (f) 
visibilities expected with an ideal tritter W' . 
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FIG. 4. Experimental reconstruction of the unitary matrix U' r of the implemented tritter at A = 785 nm compared with the theoretical one 
U 1 . (a)-(b) Real and imaginary parts of U 1 r . (c)-(d) Real and imaginary parts of U l . 



The results of the reconstruction are shown in Figs. [3]and[4] leading to the following reconstructed matrix: 

/ 0.656 ±0.001 0.5439 ±0.0006 0.5233 ± 0.0006 \ 

U r = 0.5302 ± 0.0006 (0.6135 ± 0.0006)e l ( 2 - 210±0 002 ) (0.5852 ± 0.0004)e- I ( 2 077±0 001 ) (22) 
\0.5371 ± 0.0007 (0.5725 ± 0.0003) e - l ( 2128±0 001 ) (0.6194 ± 0.0008)e l ( 2188±0 001 ). / 

We then evaluated the similarity <S' r ,m between the measured visibility matrix V m [Fig. [3(a)] and the reconstructed visibility 
matrix V r [Fig. [3(b)], leading to a value S' r . m = 0.976 ± 0.001. The correspondence between the reconstructed matrix 
U r and the ideal matrix U 1 of Eq. (O is evaluated from the similarity <S' r ,t, leading to S r ^t — 0.9595 ± 0.0003. Hence, the 
reconstructed matrix keeps a high agreeement with an ideal tritter also for 785 nm input photons. 

IV. MODELING THE EXPERIMENT 

In this section we discuss a theoretical model to take into account photon distinguishability and six photon terms in the 
theoretical prediction of the three-photon interference visibilities. Furthermore, we explicitly evalute the classical bounds for the 
Hong-Ou-Mandel visibilities by considering three independent, phase-randomized, input coherent states in the tritter. 

A. Photon distinguishability and six-photon states 

Due to the spectral properties of parametric down-conversion sources, photons belonging to different photon pairs are not 
completely indistinguishable. This effect is responsible of a reduction in the interference visibility. To include this effect in the 
theoretical predictions of the experimental visibilities, we considered as input state in the tritter the following density matrix: 

g = r\l, 1, 1)(1, 1, 1| + (1 - r)|l„, h, 1«)<1«, l b , l a \- (23) 

Here, r is the parameter which takes into account the distinguishability between the three photons, and the indexes i = a,b label 
the photon on input port 2 belongs to a different photon pair with respect to the photons on input ports 1 and 3. The parameter 
r is proportional to the overlap between the spectral functions of the three photons, and can be written as r — p 2 , where p is 
the distinguishability of two photons belonging to different photon pairs. The latter has been measured directly by performing a 
Hong-Ou-Mandel two-photon interference in a symmetric 50/50 beam-splitter, leading to p — 0.63 ± 0.03. 

Furthermore, parametric down-conversion sources in this regime present a non-negligible probability of generating three 
photon pairs. The overall state produced up to third order in the parametric gain g is given by: 

~ |0) + V2g\^ 2 ) + V3g 2 \^) + 2 5 3 |V> 6 ), (24) 

where | -0 2 ) ^ \ipi}i an d \if>e) correspond respectively to the two-photon, four-photon and six-photon contributions. In the present 
implementation, g reaches a value of g ~ 0.12. We then included in our theoretical prediction for the Hong-Ou-Mandel 
visibilities the contribution given by the six-photon term l^e)- 

B. Classical bounds for the Hong-Ou-Mandel visibilities 

We are now interested in studying the behaviour the experimentally realized tritter when it is injected by classical light (i.e. 
independent phase randomized coherent states), in order to see whether or not the obtained three photons visibilities obtained 
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by injecting a three single photons state can be ascribed to a quantum behaviour. To this end we start writing the input output 
relations as a function of the reconstructed transfer matrix W: 



b\ = YjAl^ (25) 



with i = 1,2, 3, and consider the evolution of the input state |ai)|a 2 )|a 3 ) = D(a 1 )D(a2)D(a 3 )\0), where |a^) is a generic 
coherent state and D(ai) the corresponding displacement operator. By definition: D(ai) — e ^\- a X a \ where a 2 = ^^e 101 . 
By expanding the a\ operators as a function of the bi ones, we can write the output state as: 

3 

|^ t )=IJ2e-^l a / 2 ^ T &t n '|0), (26) 

i=l rii 

where 7, = J^ij i(W r )~ 1 ]ji a j- The output probabilities related to the (300, 111, 210) output configurations can then be obtained 
by projecting \ip out ) on the respective output state, we then obtain: 

p( 3 )(lll) = e-n =1 l7,l 2 TJ| 7 .|2 ; (27) 

P( 3 )(210) = l e -^ 1 l7 1 | 2 | 7l |4| 72 |2 ; (2g) 

P( 3 )(300) = |e-n=il7 i | 2 | 7l |6. (29) 

these values correspond to the three photons perfect interference, being the three coherent states simultaneous inside the tritter. 
In order to find the corresponding dip (or peak) visibilities, we have to quantify the values of probabilities related to the delayed 
case in which one (or two) of the three coherent state is (are) delayed. If only one input state is delayed we obtain: 

P^(m u m 2 ,m 3 ) = e ~f Ll ,'^ fl U\ 2 + KTP , (30) 
mi!TO 2 !m 3 ! J -- L v ' 

z— 1 

where 

• for delayed state in port 2: 7 < = E^i^P'T 1 ]^ and %" = [(«' "T'k^, 

• for delayed state in port 3: 7- = Ylj=i ,&3 W r )~ 1 \ji a i and7 4 " = (W"') _1 ] 3i a3. 
When two over three coherent states are delayed: 

pV(m u m 2 ,m 3 ) = gg^ f[ + KT + ^TP , OD 

where 7- = [(W r ) _1 ]iiai, 7" = [(W r ) _1 ]2ia;2,7i" = [(W r ) _1 ]3ia3. Finally we perform an average over the three phases 
i = 1,2, 3, since we are dealing with phase randomized coherent states: 

/>2,7r />27r />2,7T 

r ( * ) (mi 1 m 2l ms) = ^jj y y P^{m l7 m 2 ,m 3 )de 1 d9 2 d6 3 . (32) 
From these r^ fc ^ values, we can obtain the classical visibilities: 



y = rw^ (33) 

-t OO 



where: 

• when one port is delayed: = r^ 2 ' (mi, m 2 , m 3 ) and r = r' 3 ' (mi,m 2 , 77J3) 

• when two ports are delayed: — r( 1 ^(mi, m 2 , m 3 ) and r = f^ 2 ' (mi, m 2 , 7713). 
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